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Abstract 



The B -B oscillation frequency Am^ is measured by the LHCb experiment 
using a dataset corresponding to an integrated luminosity of l.Ofb -1 of proton- 
proton collisions at yfs = 7TeV, and is found to be Arrid = 0.5156 ± 0.0051 (stat.) ± 
0.0033 (syst.) ps -1 . The measurement is based on results from analyses of the decays 
B°-> D-7T+ (D~ -> K + ir~ir~) and B° -> J/^K*° (J/0 -> K*° -> JT+7T-) 

and their charge conjugated modes. 
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1 Introduction 



The frequency Am^ of oscillations between B° mesons and B° mesons also describes the 
mass difference Am^ between the physical eigenstates in the B°-B° system, and has been 
measured at LEP II], the Tevatron |2], and the B factories (3j [4j . The current world 
average is Am^ = 0.507 ± 0.004 ps _1 15], whilst the best single measurement prior to this 
Letter is by the Belle experiment, Arrtd = 0.511 ± 0.005 (stat.) ± 0.006 (syst.) ps _1 |4J. In 
this document the convention h = c = 1 is used for all units. 

With increasing accuracy of the measurement of Am s , the counterpart of Am^ in the 
B®-B° s system j6], a more precise knowledge of Am^ becomes important, as the ratio 
Am d /Am s together with input from lattice QCD calculations [7j[8] constrains the apex 
of the CKM unitarity triangle [9,10 . Therefore, the measurement of Am^ provides an 



important test of the Standard Model 11 12 . Furthermore, Am^ is an input parameter 
in the determination of sin 2/3 at LHCb |13|. 

This Letter presents a measurement of Arrtd, using a dataset corresponding to l.Ofb -1 
of pp collisions at y/s = 7TeV, using the decay channels B° — > D~n + (D~ — > K + -k~'k~) 
and B° — > J/tpK*° (J/ip — >■ K*° — > K + n~) and their charge conjugated modes. 

For a measurement of Arrid, the flavour of the B° meson at production and decay 
must be known. The flavour at decay is determined in both decay channels from the 
charge of the final state kaon; contributions from suppressed B° — > D + n~ amplitudes 
are negligible. The determination of the flavour at production is achieved by the flavour 
tagging algorithms which are described in more detail in Sect. |4j 

The B° meson is defined as unmixed (mixed) if the production flavour is equal (not 
equal) to the flavour at decay. With this knowledge, the oscillation frequency Am^ of the 
B° meson can be determined using the time dependent mixing asymmetry 

^f(t) = ^" illZ miXe il = cos(Am,t), (1) 

'unmixed I, £ J ~r I 'mixed I, •*) 

where t is the B° decay time and -/V( un ) mixe d is the number of (un)mixed events. 



2 Experimental setup and datasets 

The LHCb detector (l4| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift-tubes placed downstream. The combined tracking 
system has a momentum resolution Ap/p that varies from 0.4% at 5 GeV to 0.6% at 
100 GeV, and an impact parameter (IP) resolution of 20 \xm for tracks with high transverse 
momentum. Charged hadrons are identified using two ring-imaging Cherenkov detectors. 
Photon, electron and hadron candidates are identified by a calorimeter system consisting of 
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scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers. The trigger consists of a hardware stage, based on 
information from the calorimeter and muon systems, followed by a software stage which 
applies a full event reconstruction. 

Events including B° — y D~n + decays are required to have tracks with high transverse 
momentum p^ to pass the hardware trigger. The software trigger requires a two-, three- or 
four-track secondary vertex with a large sum of the of the tracks, significant displacement 
from the associated primary vertex (PV), and at least one track with j?t > 1.7 GeV and 
a large impact parameter with respect to that PV, and a good track fit. A multivariate 



algorithm is used for the identification of the secondary vertices 15 . 

Events in the decay B° — y J/ip K*° are first required to pass a hardware trigger which 



selects a single muon with p T > 1.48 GeV. In the subsequent software trigger 15 , at least 



one of the final state particles is required to have pr > 0.8 GeV and a large IP with respect 
to all PVs in the event. Finally, the tracks of two or more of the final state particles are 
required to form a vertex which is significantly displaced from the PVs in the event. 
For the simulation studies, pp collisions are generated using Pythia 6.4 (16) with a spe- 



cific LHCb configuration 117] . Decays of hadronic particles are described by EvtGen 18 



in which final state radiation is generated using Photos 19 . The interaction of the 



generated particles with the detector and its response are implemented using the Geant4 



toolkit 20 as described in Ref. [21]. 



3 Selection 

The decay time t of a B° candidate is evaluated from the measured momenta and from a 
vertex fit that constrains the B° candidate to originate from the associated PV (22]| . No 
mass constraints on the intermediate resonances are applied. For the calculation of the 
invariant mass m, no mass constraints are used in the B° — y D~n + channel, while the J/ip 
mass is constrained to the world average [5] in the analysis of the decay B° — y J/ipK*°. 

All kaons, pions and muons are required to have large pt and well reconstructed tracks 
and vertices. In addition to this, particle identification is used to distinguish between pion, 
kaon and proton tracks. 

The B° — y D~n + selection requires that the D~ reconstructed mass be in a range of 
±100 MeV around the world average 15). Furthermore, the D~ decay vertex is required to 
be downstream of the PV associated to the B° candidate. 

The sum of the D~ and 7r + px must be larger than 5 GeV. The B° candidate invariant 
mass must be in the interval 5000 < m(K + 7T~7T~7T + ) < 5700 MeV. Additionally, the cosine 
of the pointing angle between the B° momentum vector and the line segment between PV 
and secondary vertex is required to be larger than 0.999. 

Candidates are classified by a boosted decision tree (BDT) |23| with the AdaBoost 



algorithm 24 . The BDT is trained with B° s — y D s 7r + candidates with no particle ID 



criteria applied to the daughter pions and kaons. The cut on the BDT classifier is optimised 
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in order to maximise the significance of the B° — y D~tt + signal. Several input variables are 
used: the IP significance, the flight distance perpendicular to the beam axis, the vertex 
quality of the B° and the D~ candidate, the angle between the B° momentum and the 
line segment between PV and B° decay vertex, the angle between the D~ momentum 
and the line segment between PV and the D~ decay vertex, the angle between the D~ 
momentum and the line segment between the B° decay vertex and D~ decay vertex, the 
IP and pt of the ir + track, and the angle between the 7r + momentum and the line segment 
between PV and B° decay vertex. Only B° candidates with a decay time t > 0.3 ps are 
accepted. 

To suppress potential background from misidentified kaons in D~ — y K~K + n~ decays, 
all D~ candidates are removed if they have a daughter pion candidate that might pass 
a loose kaon selection and are within a ±25 MeV mass window (the D~ mass resolution 
is smaller than 10 MeV) around the D~ mass when that pion is reconstructed under the 
kaon mass hypothesis. 

The B° — y J/ip K*° selection requires that the K*° candidate has a > 2 GeV and 
826 < m(K + Tt-) < 966 MeV. 

The unconstrained invariant mass must be within ±80 MeV of the J/ip mass J5]. 

B° candidates are required to have a large IP with respect to other PVs in the event 
and the B° decay vertex must be significantly separated from the PV. Additionally, B° 
candidates are required to have a reconstructed decay time t > 0.3 ps and an invariant 
mass in the range 5230 < m(J/ip K + n~) < 5330 MeV. To suppress potential background 
from misidentified B° s — y J/4"i> decays, all candidates are removed for which the K + n~ 
mass is within a ±10 MeV window around the nominal 0(1020) mass when computed 
under the kaon mass hypothesis for the pion. The resulting mass distributions for the two 
decay channels are shown in Fig. [TJ 



4 Flavour tagging 

This analysis makes use of a combination of opposite side taggers and the same side pion 
tagger to determine the flavour of the B° meson at production. The opposite side taggers, 
which use decay products of the b quark not belonging to the signal decay, are described 



in detail in Ref. 25 



The same side pion tagger uses the charge of a pion that originates from the fragmen- 
tation process of the B° meson or from decays of charged excited B mesons. Pion tagging 
candidates are required to fulfil criteria on and particle identification, as well as their 
IP significance and the difference between the B° candidate mass and the combined mass 
of the B° candidate and the pion 26 . 



Depending on the tagging decision, a mixing state q is assigned to each candidate, to 
distinguish the unmixed (q = ±1) from the mixed (q = —1). Untagged events (q = 0) 
are not used in this analysis. The tag and its predicted wrong tag probability r\ c are 
evaluated for each event using a neural network calibrated and optimized on B + — y J/ipK + , 
B ^ J/ipK*° and B° ->■ D^fi+u^ events. 
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Figure 1: Distribution of the B° candidate mass (black points) for (left) B° — » D~tt + 
and (right) B° — > J/tpK* candidates, with the results of the fits described in Sect. [6] 
superimposed. The blue dashed line shows the fit projection of the signal, the dotted 
orange line corresponds to the combinatorial background with long lifetime and the dash 
dotted red line shows the combinatorial background with short lifetime. The black solid 
line corresponds to the fit projection. 



To take into account a possible difference in the overall tagging performance between 
the calibration channels and the decay channels used in this analysis, the corrected wrong 
tag probability u assigned to each event is parametrised as a linear function of r) c (as 
described in Ref. |25| ) 

w{nc\Po,Pi) =Po+Pi(Vc~ (Vc)), (2) 

where po and p\ are free parameters in the fit for Am^ described in Sect. |6} In this way, 
uncertainties due to the overall calibration of the tagging performance are absorbed in the 
statistical uncertainty on Anid returned by the fit. 



5 Decay time resolution and acceptance 

The decay time resolution of the detector is around 0.05 ps [27]. This is small compared 
to the B° oscillation period of about 12 ps and does not have significant impact on the 
measurement of Am^. The resolution is accounted for by convolving a Gaussian function 
G(t;a t ), using a fixed width a t = 0.05 ps, with the signal probability density function 
(PDF) from Eq. ([5]). Possible systematic uncertainties introduced by the resolution are 
discussed in Sect. [3 

Trigger, reconstruction and selection criteria introduce efficiency effects that depend on 
the decay time. While these effects cancel in the asymmetry of Eq. for signal events, 
they can be important for event samples that include background. As will be shown in 
Sect, pi the only relevant background in the B° signal region is combinatorial in nature. 
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For this background the asymmetry iV (t) — N^^t^t) is expected to cancel to first order 
as q has no physical meaning. Therefore, 

, (gig) + <=1 (g) - ggUg) + <=-iW) , A 

AnixW oc (]v ^ (t) + + + oc + cos(Am d t), W 

where ^^^(t) denotes the number of unmixed or mixed signal (sig) and background 
(bkg) events. S(t) and denote the number of signal and background events as a 
function of the decay time. Thus, the shapes of S(t) and B(t) have to be known to account 
for the time dependent amplitude of the asymmetry function. 

In the analysis of B°— > J/i/jK*°, the decay time acceptance is determined from data, 
using a control sample of events that is collected without applying any of the decay time 
biasing selection criteria. The decay time acceptance is evaluated in bins of t and is 
implemented in the fit described in Sect. |6j 

In the decay B° — > D~n + there is no control dataset that can be used to measure the 
decay time acceptance. From an analysis of simulated events, it is determined that the 
decay time acceptance can be described by the empirical function 

£acc(*|oi) a 2) = arctan(ai exp(a2t)), (4) 

where the parameters a\ and a 2 are both free in the maximum likelihood fit for Am^ 
described in Sect. [6] 



6 Measurement of Am^ 

The value of Arrid is measured using a multi-dimensional extended maximum likelihood fit. 
The D~7i + data are described by a two component PDF in which one component 

describes the signal and the other describes the combinatorial background. The signal 
component consists of the sum of a Gaussian function and a Crystal Ball function 128) 
with a common mean for the mass distribution, multiplied by a function V l slg to describe 
the decay time distribution, 

Vl ig (t, q; r, Am d , u, a t , a x , a 2 ) oc 
0(t) ■ e" (1 + q(l - 2c^ c |p ,Pi))cos(Am d t)) ® G(t;a t ) ■ e^c^di, a 2 ). (5) 

Here, 0(£) is the step function, while the B° lifetime r is a free fit parameter and the 
average decay time resolution at is fixed. Other fit parameters are a\ and a 2 from the 
decay time acceptance function e acc (t|ai, a 2 ) described in Sect. [5j as well as the parameters 
Po and pi from the tagging calibration function ui{r) c \po,Pi) described in Sect. |1J Any 
B°/B° production asymmetry cancels in the mixing asymmetry function, and is neglected 
in this analysis. 

The combinatorial background component consists of an exponential PDF describing 
the mass distribution and the same decay time PDF as for the signal. For this background, 
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the decay time acceptance function is removed, Am^ is fixed to zero, and a single parameter 
is used to describe an average wrong tag probability instead of the wrong tag probability 
assigned to each event. This last parameter allows the PDF to reflect a possible asymmetry 
in the number of events tagged with q = ±1 in the background. The effective lifetime of 
the long-lived background component is allowed to vary independently in the fit. 

Possible backgrounds from misidentified or partially reconstructed decays are studied 
using mass templates determined from simulation. These are found to be negligible in the 
mass window 5200 < m(K + iT~iT~iT + ) < 5450 MeV that is used in the fit (c.f. Fig. [I]). 

In the B° — > J/ipK*° analysis, the signal mass distribution is modelled by a double 
Gaussian function with a common mean and the decay time PDF is the same as described 
in Eq. The mass distribution of the combinatorial background in B° — > J/ipK*° decays 
is also described by an exponential function. However, the decay time distribution includes 
a second component of shorter lifetime to account for prompt J/ip candidates passing the 
selection. The long-lived component is described by the same function as the combinatorial 
background in > D~it + decays, whereas the short-lived component is described by a 
simple exponential function. No other significant source of background is found. 

The fit results for B° — > D~n + and B° — > J/ipK*° decay modes are summarized in 
Sect. @ The fit yields 87 724 ± 321 signal decays for B° ->• D~t: + and 39148 ± 316 
signal decays for B°— > J/i/jK*°. The fit projections onto the decay time distributions are 
displayed in Fig. [2] and the resulting asymmetries are shown in Fig. [3j 



10' 



io 2 



10 



1 r 



.ii..! 


.... 1 . 

LHCb 




i i |. 




->t-B ->D k + 




— signal 






Ions 


b 


J = 




Vv, — combined 


: . . . 1 


i If 1 







10 15 
B° decay time t [ps] 



— , , , . , -g 

LHCb e 
+ B°^J/ V K*° Z - 

— signal 

— - short bg \ 
long bg j 

— combined 




10 15 
B° decay time t [ps] 



Figure 2: Distribution of the decay time (black points) for (left) B° — > D~tt + and (right) 
B° — > J/ipK*° candidates. The blue dashed line shows the fit projection of the signal, the 
dotted orange line corresponds to the combinatorial background with long lifetime and the 
dash dotted red line shows the combinatorial background with short lifetime (only in the 
B° — > J/ipK*° mode). The black solid line corresponds to the projection of the combined 
PDF. 
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Figure 3: Raw mixing asymmetry A m i x (black points) for (left) B° — » D~n + and (right) 
B° — > J/ipK*° candidates. The solid black line is the projection of the mixing asymmetry 
of the combined PDF. 



7 Systematic uncertainties 

As explained in Sect. [5j systematic effects due to decay time resolution are expected to 
be small. This is tested using samples of simulated events that are generated with decay 
time distributions given by the result of the fit to data and convolved with the average 
measured decay time resolution of 0.05 ps. The event samples are then fitted with the 
PDF described in Sect. [6], with the decay time resolution parameter fixed either to zero 
or to a t = O.lOps. The maximum observed bias on Am^ of 0.0002 ps -1 is assigned as 
systematic uncertainty. Systematic effects due to decay time acceptance are estimated in 
a similar study, generating samples of simulated events according to the nominal decay 
time acceptance functions described in Sect. [5} These samples are then fitted with the 
PDF described in Sect. |6j but neglecting the decay time acceptance function in the fit. 
The average observed shift of 0.0004 ps" 1 (0.0001 ps^ 1 ) in B° D~-k + (B° J/ipK*°) 
decays is taken as systematic uncertainty. 

In order to estimate systematic effects due to the parametrisation of the decay time 
PDFs for signal and background, an alternative parametrisation is derived with a data- 
driven method, using sWeights (29] from a fit to the mass distribution. The sWeighted 
decay time distributions for the signal and background components are then described 
by Gaussian kernel PDFs, which replace the exponential terms of the decay time PDF. 
This leads to a description of the data which is independent of a model for the decay time 
and its acceptance, that can be used to fit for Am^. The resulting shifts of 0.0037 ps -1 
(0.0022 ps -1 ) in the decay B° -)■ D~n + (B° -)■ J/ipK*°) are taken as the systematic 
uncertainty due to the fit model. 

Uncertainties in the geometric description of the detector lead to uncertainties in the 
measurement of flight distances and the momenta of final state particles. From alignment 
measurements on the vertex detector, the relative uncertainty on the length scale is known 
to be smaller than 0.1 %. This uncertainty translates directly into a relative systematic 



7 



uncertainty on Am d , yielding an absolute uncertainty of 0.0005 ps -1 . 

From measurements of biases in the reconstructed J/ip mass in several run periods, the 
relative uncertainty on the uncalibrated momentum scale is measured to be smaller than 
0.15%. This uncertainty, however, cancels to a large extent in the calculation of the B° 
decay time, as it affects both the reconstructed B° momentum and its reconstructed mass, 
which is dominated by the measured momenta of the final state particles. The remaining 
systematic uncertainty on the decay time is found to be an order of magnitude smaller 
than that due to the length scale and is neglected. 

A summary of the systematic uncertainties can be found in Table [I] The systematic 
uncertainty on the combined Arrid result is calculated using a weighted average of the 
combined uncorrelated uncertainties in both channels. The uncertainty on the length scale 
is fully correlated across the channels and therefore added after the combination. 

Table 1: Systematic uncertainties on Am^ in ps _1 





B ^ J/ipK* 


B ^ D--K+ 


Acceptance 


0.0001 


0.0004 


Decay time resolution 


0.0002 


0.0002 


Fit model 


0.0022 


0.0037 


Total uncorrelated 


0.0022 


0.0037 


Length scale 


0.0005 


0.0005 


Total including correlated 


0.0023 


0.0037 



8 Conclusion 

The B°-B° oscillation frequency Am d has been measured using samples of B° — > D~tt + 
and B°— > J/vJjK*° events collected in 1.0 fb _1 of pp collisions at y/s = 7TeV and is found 
to be 

Am d (B° IT 7T+) = 0.5178 ± 0.0061 (stat.) ± 0.0037 (syst.) ps" 1 and 
Am d (B ^ J/ipK* ) = 0.5096 ±0.0114 (stat.) ± 0.0022 (syst.)ps^ 1 . 

The combined value for Am d is calculated as the weighted average of the individual results 
taking correlated systematic uncertainties into account 

Am d = 0.5156 ± 0.0051 (stat.) ± 0.0033 (syst.) ps -1 . 

It is currently the most precise measurement of this parameter. The relative uncertainty 
on Am d is 1.2%, where it is around 0.6% for Am s (6J. Thus, the uncertainty on the ratio 
Am d /Am s is dominated by Am d . As the systematic uncertainties in the Am d and Am s 
measurements are small, the error on the ratio can be further improved with more data. 
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